Magnetic cryocooling with Gd^"*" centers in a light and compact framework 
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The magnetocaloric effect of gadolinium formate, Gd(00CH)3, is experimentally de- 
termined down to sub-Kelvin temperatures by direct and indirect methods. This 3D 
metal-organic framework material is characterized by a relatively compact crystal 
lattice of weakly interacting Gd^^ spin centers interconnected via light formate lig- 
ands, overall providing a remarkably large magnetic:non-magnetic elemental weight 
ratio. The resulting volumetric magnetic entropy change is decidedly superior in 
Gd(00CH)3 than in the best known magnetic refrigerant materials for liquid-helium 
temperatures and low-moderate applied fields. 
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Recent years have witnessed a terrific increase in the number of molecule-based materials 
proposed as magnetic refrigerants for liquid-helium temperatures.-"— Refrigeration proceeds 
adiabatically via the magnetocaloric effect (MCE), which describes the changes of magnetic 
entropy {ASm) and adiabatic temperature (AT^rf), following a change in the applied mag- 
netic field (AS). As in the first paramagnetic salt that permitted sub-Kelvin temperatures 
to be reached in 1933/^ gadolinium is often present because its orbital angular momentum 
is zero and it has the largest entropy per single ion.^ How to spatially assemble the Gd^^ 
spin centers is vital for designing the ideal magnetic refrigerant On the one hand, the 
magnetic density should be maximized, for instance, by limiting the amount of non-magnetic 
elements which act passively in the physical process. On the other hand, the magnetic or- 
dering for B = should not develop, causing the decrease of MCE, above the target working 
temperature of the refrigerant. Therefore a compromise becomes necessary, especially for 
reaching very low temperatures. 

This letter focuses on gadolinium formate, whose chemical formula reads Gd(00CH)3, 
which belongs to the class of metal-organic framework (MOP) materials. Other three- 
dimensional MOP materials have recently attracted the interest for their cooling properties, 
combined with their synthetic variety and intrinsic robustness.—"— We shall see below that, 
while presenting a sub-Kelvin ordering temperature, Gd(00CH)3 is also characterized by 
a relatively high packing density of Gd^+ ions, linked only to light formate ligands. The 
resulting MCE, that we estimate by direct and indirect methods, sets Gd(00CH)3 in an 
enviable position within this research area. 

Single-crystal structure determination of Gd(00CH)3 completes the original results de- 
rived from powder diffraction.— No previous magnetic measurements on Gd(00CH)3 
are reported in the literature, except for initial Mossbauer experiments.— Magnetization 
measurements down to 2 K and heat capacity measurements using the relaxation method 
down to ~ 0.35 K were carried out on powder samples by means of commercial setups for 
< i? < 5 T and < i? < 7 T, respectively. Direct measurements of the MCE were 
performed on a powder sample using a dedicated thermal sensor, installed in the same setup 
employed for the heat capacity experiments. 

Pigure 1 shows the measured molar magnetization M for temperatures within 2 — 10 K. 
The magnetization saturates to the expected value of 7 /ib for a Gd'^'^ spin moment, according 
to which s = 7/2 and g = 2. The M[T) curves can be well described by a Brillouin function 
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- see, e.g., the dashed hne in Fig. 1 for an ideal paramagnet at T = 2 K. Deviations of the 
experimental data from the paramagnetic behavior are barely noticeable only for the lowest 
temperatures, and can be ascribed to the presence of a weak antiferromagnetic interaction. 
This is corroborated by the T-dependence of the magnetic susceptibility x- As shown by the 
solid line in the inset of Fig. 1, the susceptibility data can be fitted above 2 K to a Curie- 
Weiss law X = g^fJ'B^i^ + l)/[3A;B(r — 0)], obtaining a negative, though small, 6 = —0.3 K, 
which suggests that the Gd'^+ moments are weakly antiferromagnetically correlated in the 
paramagnetic phase. 

The top panel of Figure 2 shows the measured low-temperature heat capacity C, nor- 
malized to the gas constant R, as a function of temperature for several applied fields. A 
sharp lambda-like peak can be observed in the zero-field data for Tqi — 0.8 K, denoting 
the presence of a phase transition, which is accompanied by a smooth and tiny feature at 
Tc2 ~ 0.4 K. The magnetic origin of both anomalies is proved by the fact that external 
applied fields quickly and fully suppress them.— In agreement with M(T, B), the analysis of 
the field-dependent C reveals that magnetic interactions between the Gd^+ spin centers are 
relatively weak, since applied fields i? > 1 T are sufficient for fully decoupling all spins. As 
indeed shown in Fig. 2, the calculated Schottky contributions (solid lines) for the field-split 
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FIG. 1. Field-dependence of the experimental molar magnetization M for temperatures ranging 
from 2 to 10 K, with a 1 K step between adjacent isothermal curves. Dashed line is the calculated 
M of a paramagnet for s = 7/2, g = 2 and T = 2 K. Inset: temperature-dependence of the inverse 
of the experimental molar susceptibility x collected for B = 0.1 T and Curie- Weiss fit (solid line). 
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FIG. 2. Top: temperature-dependence of the heat capacity C normahzed to the gas constant 
R collected for B = 0,1,3, and 7 T, as labeled. Solid thick lines are the calculated Schottky 
contributions for the corresponding B, and dashed line is the fitted lattice contribution. Bottom: 
T-dependence of the experimental magnetic entropy Sm normalized to the gas constant R for 
several B, as obtained from the magnetic contribution Cm to the total heat capacity. Highlighted 
examples of magnetic entropy change ASm between states A -H- D and adiabatic temperature 
changes ATad between states A o C and B o D. 

levels of the non-interacting s = 7/2 multiplet nicely account for the magnetic contribution 
Cm to the experimental heat capacity. For T > 7 K, a large field-independent contribution 
appears, which can be attributed to the lattice phonon modes of the crystal. The dashed 
line in the top panel of Fig. 2 represents a fit to this contribution, with the well-known Debye 
function yielding a value of 0d = 168 K for the Debye temperature, which is remarkably 
large among molecular— and MOF^^ materials, denoting a relatively stiff lattice. A larger 
Bd implies a correspondingly lower lattice entropy in the low-T region, ultimately favoring 
the MCE. From the experimental heat capacity the temperature dependence of the magnetic 

4 



T 1 ' 1 ' 1 ' 1 ' 1 ' T 




5 10 15 20 25 30 



r/K 



FIG. 3. Top: temperature-dependence of the magnetic entropy change ASm, as obtamed from 
magnetization and heat capacity data (Figs. 2 and 3, resp.) for the indicated apphed-field changes 
AB. Vertical axis reports units in J kg^^ K^^ (left) and volumetric mJ cm^^ (right). Bottom: 
T-dependence of the adiabatic temperature change ATad, as obtained from heat capacity data for 
the indicated AB. 

entropy Sm(T) is derived by integration, i.e., 

SUT) = l^^^dT, (1) 

where Cm is obtained by subtracting the lattice contribution to the total C measured. The 
so-obtained Sm{T) is shown in the bottom panel of Fig. 2 for the corresponding applied 
fields. For B = 0, the lack of experimental Cm for T ^ 0.3 K has been taken into account 
by matching the limiting Sm at high T with the value obtained from the in-field data. One 
can notice that there is a full entropy content of i?ln(8) per mole Gd^"*" involved, as expected 
from R\n{2s + 1) and s = 7/2. 

Next, we indirectly evaluate the MCE of Gd(00CH)3 from the experimental data pre- 
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FIG. 4. Time evolution of the applied field B, experimental temperature T and deduced adiabatic 
temperature Tad, as labeled, during a magnetization (top) and a demagnetization (bottom) process, 
both starting from Bi and Tj = 0.95 K. 

sented so far. From the bottom panel of Fig. 2, we straightforwardly obtain the mag- 
netic entropy changes AS'm(T, Ai?) for different apphed field changes /S.B = Bj — Bi. 
The so-obtained results are depicted in Figure 3. A similar set of data can also be de- 
rived from an isothermal process of magnetization by employing the Maxwell relation, i.e., 
ASm{T, AB) = /|/[9M(T, B)/dT]BdB. From the experimental M(T, B) data in Fig. 1, we 
then obtain curves that rather beautifully corroborate the corresponding results previously 
derived from heat capacity - see top panel of Fig. 3. Furthermore, to a cooling process 
under adiabatic conditions, one naturally associates a temperature change whose estimate is 
made feasible by knowing C and thus Sm- The bottom panel of Fig. 3 shows ATad{T, AB), 
where T denotes the final temperature of the adiabatic cooling, e.g., for going from states 
C(T = 3.4 K, 5 = 1 T) to A(T = 0.95 K, E = 0) in Fig. 2. 



A far more elegant and reliable method for determining the MCE is by directly mea- 
suring ATadiT, AB) under quasi-adiabatic conditions. The procedure comprises a full 
magnetization-demagnetization cycle, during which the experimental T and B are contin- 
uously recorded. In a half cycle, starting with the sample at an initial Tj, we magnetize 
(demagnetize) by gradually increasing (decreasing) the applied field from Bi to Bj and let 
the sample relax to the final Tf. In order to compute the temperature evolution for an 
ideal adiabatic process, one requires a precise knowledge of the heat that unavoidably is 
absorbed from (released to) the thermal bath during the direct measurement. For this pur- 
pose, the thermal conductance k of the wires holding our sensor was previously determined 
as a function of T, using a standard copper piece as the sample. The non-adiabaticity in- 
duces a variation of the entropy AS* = S{t) — S{tQ) in a time interval t — Iq, which can 
be expressed as AS* = //^ A;(Tj — T)/Tdt at every time instant. Besides, from Eq. [1] we 
also have AS" = J^^^C{T, B)/TdT, where the adiabatic temperature Tad, viz., the temper- 
ature if the sample would have been kept thermally isolated, is the only unknown and can 
therefore be deduced numerically. In our treatment, we safely disregard the entropy contri- 
bution due to the heat transferred from the sample holder to the refrigerant material, i.e., 

rp 

ASsh = St! ^sh/TdT, since the heat capacity of the sample holder Csh is negligible with 
respect to that of Gd(00CH)3 below liquid-helium temperature. 

Figure 4 shows the time evolution of B,T and Tad for a full magnetization-demagnetization 
cycle, starting at Tj = 0.95 K and for a field change AB = B f — Bi = {1 — 0) T or {0 — 1) T, 
depending on whether we deal with the magnetization or demagnetization process, respec- 
tively. We note that the exact same conditions are highlighted in Figure 2: process A — C 
for the magnetization and process D — B for the demagnetization. In the top panel of 
Fig. 4, we observe T to increase while we magnetize to 1 T. Here Tad increases more than T 
because the thermal losses to the bath are compensated to obtain Tad- Upon reaching Bf, 
T decays back towards Tj = 0.95 K but Tad = 3.5 K is constant, since it corresponds to an 
adiabatic process at constant B. In the bottom panel, T decreases below Tj, while we de- 
magnetize to zero field, whereupon T gradually relaxes back to equilibrium, while constant 
Tad = 0.45 K. Remarkably, the final adiabatic temperatures of 3.5 K and 0.45 K obtained 
after sweeping the field up and down, respectively, corroborate the results independently 
inferred by an indirect method - see states C and B, respectively, in Figure 2. 

The MCE of Gd(00CH)3 is exceptionally large, especially in comparison with other 
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molecule-based magnetic refrigerants. A fine example is the recently studied dimeric light 
molecule [{Gd(OAc)3(H20)2}2] ■4H2O, hereafter shortened as {Gd2}, whose —ASm reaches 
a value as large as 40.6 J kg^^ K^^ for AB = (7 — 0) T and T ~ 1.8 K.^ For widespread 
applications, the interest is chiefly restricted to applied fields which can be produced with 
permanent magnets, viz., in the range 1 — 2 T. In this regard {Gd2} could be appealing 
because a weak, ferromagnetic intracluster exchange interaction enhances the MCE for low 
applied fields, yielding an outstanding —ASm = 27.0 J kg~^ for AB = (1 — 0) T and 
T ~ 0.8 K. However, these values are not as large as the ones we here report for Gd(00CH)3, 
e.g., -ASm = 56.0 J kg^i K^^ and 31.2 J kg-^ K'^ for A5 = (7 - 0) T and (1 - 0) T, 
respectively, at similar corresponding temperatures such as seen in Figure 3. 

Gadolinium gallium garnet (GGG) is the reference magnetic refrigerant material for the 
liquid-helium temperature region.— Indeed, its functionality is commercially exploited in 
spite of a relatively modest maximum —ASm = 20.5 J kg~^ for AB = (2 — 0) T. This 
apparent contradiction is resolved by measuring the entropy change in terms of equivalent 
volumetric units, which take into consideration the GGG mass density p = 7.08 g cm~'^. By 
so-doing, GGG achieves a record value —pASm — 145 mJ cm~'^ for the same applied 
field change of 2 T. Although these units are not often used to characterize the MCE, they 
are better suited for assessing the implementation of the refrigerant material in a designed 
apparatus.— On this point, one could correctly argue that the MCE of molecule-based 
refrigerant materials is disfavored by their typically low p. For instance, p = 2.04 g cm~'^ for 
the aforementioned {Gd2}, which results in —pASm — 56 mJ cm~'^ and 73 mJ cm~^ 
for AB = (1 — 0) T and (3 — 0) T, respectively,- i.e., definitely much lower than GGG. 

The mass density p = 3.86 g cm~^ of Gd(00CH)3 is very large among molecule-based 
materials, though yet smaller than that of GGG. In Gd(00CH)3, the Gd'^^ centers are 
interconnected only by short and extremely light CHOO~ ligands. Ultimately, this enhances 
the MCE favored by a larger weight of magnetic elements with respect to non-magnetic 
ones, which act passively. As a matter of fact, the mass density of these two materials is 
effectively counterbalanced by the magnetic:non-magnetic weight ratio, which amounts to 
0.54 in Gd(00CII)3 and to a lower 0.47 in GGG.^^ For comparison, this number further 
reduces to 0.39 in the case of {Gd2}. Overall, adopting the proper units, the MCE of 
Gd(00CH)3 is characterized by maxima —pASm — 120 mJ cm~'^ and 189 mJ cm~^ 
for AB = (1 — 0) T and (3 — 0) T, respectively, as can be seen in Figure 3. These values 
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compares favorably with the ones obtained from GGG. 

Concluding, we experimentally determine the magnetocaloric effect of the Gd(00CH)3 
metal-organic framework material. Under quasi-adiabatic conditions, sub-Kelvin direct mea- 
surements of the temperature change corroborate the results inferred from indirect methods. 
By comparing the MCE per volume of other known materials, such as GGG, we demonstrate 
that gadolinium formate could serve as an excellent magnetic refrigerant for liquid helium 
temperatures. Our observations are interpreted as the result of a light and compact struc- 
tural framework promoting very weak magnetic correlations between the Gd'^"'" spin centers. 
Finally, we foresee that synthetic and technological strategies, already developed for the sur- 
face deposition of MOF materials, could ultimately facilitate the integration and exploitation 
of Gd(00CH)3 within molecule-based microdevices for on-chip local refrigeration.— 
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For Gd(00CH)3, the relative molecular mass and number of Gd^"*" ions per formula unit 
are m = 292.30 g mol~^ and n = 1, respectively. For GGG, m = 1012.37 g mol~^ and 
n = 3. For {Gd2}, m = 812.89 g mol~^ and n = 2. Therefore, nAr/m = 0.54, 0.47 and 
0.39 for Gd(00CH)3, GGG and {Gda}, respectively, where Ar = 157.25 g mor^ is the 
gadolinium relative atomic mass. 
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